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Abstract 

Results of a systematic investigation of tlie vapour-liquid equilibria of 38 indi- 
vidual two-centre Lennard-Jones plus axial pointdipole model fluids (2CLJD) are 
reported over a range of reduced dipolar momentum < fi*'^ < 20 and of reduced 
elongation < L* < 1.0. Temperatures investigated are from about 55 % to about 
95 % of the critical temperature of each fluid. The A^T-I-Test Particle Method is 
used for the generation of vapour pressures, saturated densities, and saturated en- 
thalpies. For the lowest temperatures, these data are calculated with highly accurate 
chemical potentials obtained from the gradual insertion method. Critical temper- 
atures T* and densities p* are obtained from Guggenheim's equations. Empirical 
correlations for critical data T* and as well as for saturated densities p'* , p"* , 
and vapour pressures p* are developed as global functions of the model parameters. 
They describe the simulation data generally within their statistical uncertainties. 
Critical pressures and acentric factors of the 2CLJD fluid can be calculated from 
these correlations. The present results are a sound basis for adjustments of the 
model parameters p*'^, L*, a, and e to experimental VLE data of real fluids. 
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1 Introduction 



Knowledge of vapor-liquid equilibria (VLE) is important in many problems in engineering 
and natural sciences. Among the different ways to model vapor-liquid equilibria, molec- 
ular simulation has the highest potential to yield significant improvements compared to 
existing models, especially in terms of predictive power. However, further efforts are 
needed until molecular simulation based models and tools will be sufficiently developed 
so that their advantages can help process engineers in their real world tasks. 

One of the main problems to be overcome is the lack of intermolecular interaction 
models describing vapor-liquid equilibria of real fluids with technically relevant accuracy. 
As in most applications mixtures are of interest, it is necessary to model pure fluids with 
a family of compatible interaction models, allowing their application for mixtures using 
simple combining rules. 

Various Stockmayer fluids, i.e. one-centre Lennard- Jones plus pointdipole (ICLJD) 
fluids, and two-centre Lennard- Jones plus pointdipole (2CLJD) fluids have been studied 
by other authors. These model fluids have been applied to a number of real dipolar fluids. 
Stockmayer parameters are presented by van Leeuwen [1] for various real substances and 
by Gao et al. [2J for alternative refrigerants. 2CLJD parameters for various refrigerants 
are reported by Vega et al. [ISj, Kohler and van Nhu ^, Kriebel et al. |5j and by Lfsal et 
al. [6]. Model parameters for the dipolar Kihara potential, which is similar to the 2CLJD 
potential model, are given by Lago et al. [7j for organic solvents. 

The search for an appropriate interaction model for a given fluid is usually a time 
consuming process. In general, the focus lies on one pure fluid, where the optimization of 
the potential model is done by a number of simulations with subsequent variation of the 
model parameters, cf., e.g., van Leeuwen et al. [8] for methanol. 
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In a previous work we followed a new route to develop quantitative interaction models 
that allows fast adjustments of model parameters to experimental data for a given class 
of pure fluids. For the two-centre Lennard- Jones plus pointquadrupole (2CLJQ) fluid 
systematic studies of the vapour-liquid equilibria were carried out. The results were 
correlated as a function of the model parameters of the 2CLJQ fluid [9J. Using these 
correlations, it was possible to determine the model parameters for a large number of real 
quadrupolar fluids ^U\. It has been shown, that these molecular models can successfully 
be applied for the description of vapour-liquid equilibria of binary and multicomponent 
mixtures [TT] . 

Encouraged by these favorable results, in this work an analogous investigation is car- 
ried out for the important class of dipolar fluids. Only a few studies of VLE of the 2CLJD 
fluid are available in the literature [HI 1121 1131 HI] , which however cover only a compara- 
tively narrow range of the model parameters. Furthermore, VLE of some submodels of 
the 2CLJD fluid have previously been investigated: Galassi and Tildesley [I5j, Kriebel et 
al. [IE], and Kronome et al. [IT] the 2CLJ, van Leeuwen [I], van Leeuwen et al. [T8] . 
Gao et al. [2], Garzon et al. [TH], and Smit et al. ^2^ the Stockmayer fluid (ICLJD). 
The simulation techniques applied were Gibbs-Duhem integration or the Gibbs Ensemble 
Monte Carlo method. The idea followed in this work is to study the VLE of the 2CLJD 
model fluid systematically and in detail over the whole relevant range of model param- 
eters. Using reduced coordinates, for the symmetric 2CLJD fluid, only two parameters 
have to be varied: the dimensionless LJ centre-centre distance L* and the dimensionless 
squared dipolar momentum /i*^. The parameter space of interest can be covered with ac- 
ceptable accuracy by studying 38 individual 2CLJD fluids with different values for L* and 
/i*^. The simulation results are correlated empirically in order to be able to interpolate 
between the discrete pairs of L* and /x*^. 



3 



Given the results from the present study, it is straightforward to adjust the molecular 
interaction parameters of the 2CLJD fluid to experimental VLE data of real dipolar 
fluids. Properties like the critical values of temperature and pressure, the acentric factor 
or the saturated liquid density and vapour pressure are available as functions of the 
molecular interaction parameters. Hence, the development of the molecular interaction 
model for a given substance is not more difficult than the adjustment of parameters of 
phenomenological thermodynamic models. 

This work covers the basic molecular simulations for 38 individual 2CLJD fluids and 
the development of the correlations together with a discussion of these results. The 
application to real fluids will be presented separately. 

For the calculation of vapour-liquid equilibria the A^T+Test Particle Method [21] was 
chosen due to favourable experience with that method especially concerning accuracy. 
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2 Investigated model class 



In this work pure two-centre Lennard- Jones plus axial pointdipole (2CLJD) fluids are 
studied. The 2CLJD potential model is composed of two identical Lennard- Jones sites a 
distance L apart (2CLJ), forming the molecular axis, plus a pointdipole (D) n placed in 
the geometric centre of the molecule along the molecular axis. The pair potential ii2CLjD 
and further technical details are described in the Appendix A. 

The parameters a and e of the 2CLJD pair potential were used for the reduction of 
the thermodynamic properties as well as the model parameters L and jj?-. T* = Tk/e, 
p* = pcrVe, P* = h* = h/e, L* = L/a, fi*^ = /^V (47reo • ea^), At* = At^m/e/a. 

The reduced parameters L* and fi*^ were varied in this investigation: L* = 0; 0.2; 0.4; 
0.505; 0.6; 0.8; 1.0 and //*^ — 0; 3; 6; 9; 12. All combinations of these values lead to a 
set of 35 model fluids. Additionally to these, the fluids with L* — 0.0; 0.2 and /j,*^ — 16 
as well as L* = 0.0 and /i*^ = 20 were investigated, as very strong dipolar momenta are 
only reahstic for molecules with small elongation L*. So, altogether 38 model fluids are 
considered. 
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3 Molecular simulation method for VLE data 



For all 38 model fluids the NpT+Test Particle Method (A^T+TP Method) proposed 
by Moller and Fischer [HI 122 EH] was applied to obtain the VLE data. The NpT+TF 
Method performs separate NpT simulations in the liquid and the vapour phase and uses 
information on the chemical potentials to calculate the VLE. Vapour pressures, saturated 
densities and residual enthalpies h^'^^* {T*,p*) = h* (T*,p*) — h^^* (T*) in equilibrium are 
evaluated. 

Test particle insertion for the calculation of the residual chemical potential fi^^^* in 
the NpT ensemble is based in general on Widom's method [2^, which, however, yields 
with large uncertainties at dense low-temperature state points. This effect is partic- 
ularly important when highly dipolar fluids with large elongations are considered. As the 
uncertainty of f/'^^* directly influences the uncertainty of the vapour pressure, Widom's 
method would spoil any accurate calculation of vapour pressures at these state points. 
More precise values for the residual chemical potential can be obtained from the gradual 
insertion method flEl flEl fI7\ . This method was used here in order to increase the accu- 
racy of vapour pressure calculations at lowest temperatures 0.55 • T* of the systems with 
L* G [0; 1.0] and /i*^ g [0; 12], except for L* = with /i*^ = 12. For the systems L* = 
with = 12, L* = or L* = 0.2 with /i*^ = 16, and L* = with /i*^ = 20 Widom's 
method was applied with an increased number of test particles and extended simulation 
run lengths at state points 0.55 ■ T*. It must be pointed out, that due to these difficulties 
at low temperature state points, VLE simulations of dipolar fluids are performed typically 
down to about 0.70 ■ T*. Below, special techniques are required for good results. 

Configuration space sampling was done by A^T-molecular dynamics simulations with 
= 864 particles for both liquid and vapour simulations. The dimensionless integration 
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time step was set to At* = 0.0015. Starting from a face-centred lattice arrangement every 
simulation run was given 10, 000 integration time steps to equilibrate. Data production 
was performed over n = 100, 000 (n = 200, 000 for the three systems mentioned above) 
integration time steps. At each production time step 2N {ION for the three systems 
mentioned above) test particles were inserted in the liquid phase, and test particles 
were inserted in the vapour phase, in order to calculate the chemical potentials. The di- 
mensionless dynamical parameter of A^T-MD-simulations ascribed to the box membrane 
according to Andersen's algorithm [28] was set to 2 ■ 10"'' for liquid simulations and to 
10~^ for vapour simulations. The high value of A^ allowed simulations in the vicinity of 
the critical point, and the high value of n was used in order to obtain small statistical 
uncertainties. 

The gradual insertion method for the calculation of residual chemical potentials fx^^^* is 
an expanded ensemble method based on the Monte Carlo technique. The gradual insertion 
method as proposed in p6] extended to NpT ensembles [27] was used in this work. In 
comparison to Widom's test particle method, where real particles are inserted in the 
fluid, in gradual insertion one fluctuating particle is introduced, that undergoes changes 
in a predefined set of discrete states of coupling with all other real particles of the fluid. 
Preferential sampling is done in the vicinity of the fluctuating particle. This concept leads 
to considerably improved accuracy of the residual chemical potential. Gradual insertion 
simulations were performed with A^ = 864 particles in the liquid phase at 0.55 ■ T*. 
Starting from a face-centred lattice arrangement every simulation run was given 5, 000 
Monte Carlo loops to equilibrate. Data production was performed over n = 100, 000 
Monte Carlo loops. One Monte Carlo loop is deflned here as A^ trial translations, (2/3) A^ 
trial rotations, and one trial volume change. Further simulation parameters for runs with 
gradual insertion were taken from Vrabec et al. fI7\ . 
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For the 2CLJD fluid a hybrid equation of state (EOS) for the Helmholtz energy -F2CLJD 
is available [3, [291 ED] , which is constructed as sum of the Helmholtz energy -F2CL j for the 
Lennard- Jones part of the potential model and the Helmholtz energy Fd for the dipolar 
contribution. In planning the present investigation this 2CLJD EOS was used to estimate 
the critical temperatures T* (/i*^, L*) for the systems with L* G [0; 0.8] and /i*^ G [0; 12]. 

For all systems VLE data were calculated for temperatures from about 55 % to about 
95 % of the critical temperature. Liquid simulations were performed for the whole tem- 
perature range, whereas vapour simulations were only performed at temperatures above 
0.80 ■ T^*estd- Below this temperature the second virial coefficient is sufficient for the VLE 
calculations. Liquid simulations took about eleven hours, vapour simulations about three 
hours CPU time on a modern workstation (e.g., Compaq AlphaStation XPIOOO). Liquid 
simulations with gradual insertion took an average of about 200 hours CPU time. 
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4 Simulation results of VLE data and critical data 



Table [H reports an extract of the VLE data obtained in this work. All 38 fluids are 
covered, but not all state points for which simulations were carried out are included in 
Table [T] for the sake of brevity. For each state point, the vapour pressure p* , the saturated 
liquid density p'*, the saturated vapour density p"*, the residual saturated liquid enthalpy 
h'^'^^* , and the residual saturated vapour enthalpy h"^'^^* are reported. The temperatures 
are about T* ^ 0.55 -T^*, 0.80 -T^.*, and 0.95 -T^.*. Statistical uncertainties were determined 
with the method of Fincham et al. [31j and the error propagation law. The full data is 



available at http://www.itt.uni-stuttgart.de/molsim.html 



Figs. [T] - m illustrate for /i*^ = 3 and p*^ = 12 the strong influence of both the elonga- 
tion and the dipolar momentum on the 2CLJD VLE data. Increasing the elongation or 
increasing the dipolar momentum strongly influences the shape of the density coexistence 
curve and the slope of the vapour pressure curve. At low temperatures, the high uncer- 
tainties of the vapour pressures of the systems with large elongations and strong dipolar 
momentum, cf. Fig. IU are due to the uncertain values of the chemical potentials obtained 
by Widom's test particle insertion in the liquid phases. Figs. |3] and H] illustrate, at lowest 
temperatures, the considerable decrease of the uncertainties of the vapour pressure when 
the chemical potential in the liquid phase is obtained with much lower uncertainties by 
gradual insertion. 

The comparison of the 2CLJD VLE data from the present work to the 2CLJD EOS 
reveals systematic deviations. This can best be seen in Figs. [2] and HI For models with 
small dipolar momenta and for models with elongations near L* = 0.505 the 2CLJD VLE 
data from the simulations and that from the EOS agree well. This is due to the fact, 
that the dipolar contribution to the EOS is based on data of a 2CLJD model fluid with 
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L* = 0.505. However, for other state points the EOS sometimes deviates considerably 
from the simulation data. In most cases the EOS underestimates the saturated liquid 
densities and overestimates the saturated vapour densities. The vapor pressures from the 
2CLJD EOS generally agree well with simulation data from this work, except for model 
fluids with high dipolar momenta and large elongations, where the EOS overestimates the 
vapour pressures (cf. Fig. HJ the EOS is not valid beyond L* = 0.8). 

Critical data were determined here with the method of Lotfi et al. [35], who, by simple 
means, found reliable critical data for the ICLJ model fluid. It is known that the density- 

1 /3 

temperature dependence near the critical point is well described by p* ~ (T* —T*)' , as 
given by Guggenheim [331 [3l]. Eqs. ([T]) and ([2]), suggested by Lotfi et al. |32j, were used 
for the correlation of the saturated densities from simulation 

p'* = pi + Ci ■ (t; - T*)i/3 + c'^ ■ (t; - t*) + c'^ ■ {t* - t*)^/\ (i) 
p"* = pi-Ci- (t; - T*y/^ + c'i ■ (t; - t*) + c'l ■ {t* - T*fi^. (2) 

The simultaneous fit of saturated liquid and saturated vapour densities yields not only the 
coefficients Ci, C2, C3, , C3, but also the critical data p*, T*. The critical temperatures 
and densities for the 2CLJD model fluids are listed in Table [21 In most cases critical 
temperatures obtained from the simulation data are lower than those estimated by the 
2CLJD EOS. Table El also contains the critical compressibility factor = p*/ {p*T*), 
which is a non-reduced property of the 2CLJD fluids. It is therefore of particular interest 
for comparisons to real dipolar fluids. 

The uncertainties of T* and p* are similar as those found in an analogous investigation 
of the two-centre Lennard- Jones plus pointquadrupole fluid [H]. They are estimated to 
be ciT*) ^ 0.005 and cr(p*) ~ 0.0005, and it is concluded, that the critical temperatures 
calculated by this method are certain up to the second, the critical densities up to the 
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third digit after the decimal point. 



5 Global correlation of VLE data 

In order to obtain VLE data for the whole range of /x*^, L* and T* the molecular simulation 
data from this work were globally correlated. The critical data T*{fi*'^, L*), pKfi*'^, L*), the 
saturated hquid density p'*{fi*^, L*,T*) and the vapour pressure 
pl{fi*^,L*,T*) are the key VLE data for an adjustment to real fluids. The adequate 
shape of the temperature-density coexistence curve was achieved by simultaneously cor- 
relating the functions p'* (/i*^, L*,T*) and p"* {p*^, L*,T*). It was not in the scope of the 
present investigation to construct a new 2CLJD EOS. The correlation developed here is 
not designed to compete with an EOS. It shall merely be a working tool for a restricted 
field of application, namely the adjustment of model parameters to data of real fluids as it 
has been suggested by Vrabec et al. [10]. The vapour pressure correlation is also used to 
verify the thermodynamic consistency of the VLE data from simulations by the means of 
the Clausius-Clapeyron equation. Moreover, the correlations are useful for comparisons 
with results of other investigators. Details of the correlation method are described in the 
Appendix B which also contains the resulting correlation functions. 

5.1 Critical properties 

The correlation functions T*{p*'^, L*) and p^p*"^, L*) were assumed to be linear combina- 
tions of elementary functions, cf. Eq. in the Appendix B. The elementary functions 
and their coefficients are given in Table [3] in the Appendix B. The quality of the corre- 
lations can be studied in Fig. Most relative deviations of the critical temperatures 
are within 1 %. The critical densities are represented with roughly the same quality. It 
should be mentioned that possible systematic errors may be introduced to the critical 
data by the choice of the exponent 1/3 in the second term of Eqs. ([T]) and ([2]). 
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5.2 Saturated densities, vapour pressure 

The temperature-saturated density correlations are based on Eqs. ([I]) and ([2]), which 
have five adjustable parameters Ci, Cg, Cg, C2, and C3. Another three parameters are 
introduced by the correlation of the vapour pressure, cf. Eq. ([7j) in the Appendix B. Again 
the functions describing the dependency of the correlation parameters on /i*^ and L* were 
assumed to be linear combinations of elementary functions. The elementary functions 
and their coefficients are given in Table [3] in the Appendix B. 

A comparison between the correlations and the simulation data can be seen in Figs. 
[U - m A more detailed comparison is given in Fig. [HI For some of the 2CLJD fiuids 
investigated here, the relative deviations of the simulation results from the correlations 
are shown. Typically, the saturated liquid density correlation has the largest deviations 
for the highest temperature of each model, which is due to the large uncertainties in the 
near critical region. The correlation shows relative deviations in the range of 0.4 % for 
mid temperatures which are most important for adjustments to real fiuids. 

Relative deviations of the saturated vapour densities are not illustrated here. The 
vapour density correlation should not be used below 0.60-Tc(/i*^, L*). At low temperatures 
the correlation is not useful as it does not capture the limiting case of the ideal gas, which 
is independent of the parameters /x*^ and L*. 

In most cases, the vapour pressure correlation represents the simulation data within 
their uncertainties. It has to be mentioned, that, except for simulations with gradual in- 
sertion, simulations at low temperature state points yield vapour pressures with increased 
uncertainties due to the uncertain values of the chemical potential obtained by Widom's 
test particle insertion in dense liquid phases. 

By extrapolating the vapour pressure correlation slightly to the critical point, the 
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acentric factor [33] 

-(/^ '^) = -l°gio 1 (3) 

can be calculated from the correlations discussed before. Table [2] contains the critical 
pressures and the acentric factors for the 2CLJD model fluids calculated on the basis of 
these correlations. 
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6 Discussion 



6.1 Comparison to results of other authors 

The results from the present study are compared here to the simulation results of other 
authors. Fig. [7] presents results for the relative deviations for the saturated liquid density 
(Pothcr -Pcorr)/Pcorr and the vapour pressure (p^^^ther -Pa,corr)/K,corr- The saturated liquid 
densities of other authors agree in almost all cases within the combined uncertainties 
of our correlation and of the simulation results, cf. Fig. [71 top. Near critical points 
of course show larger deviations. Molecular simulations with considerably lower particle 
numbers yield slightly lower densities, thus causing systematic negative deviations, as can 
be observed, for example, for data from Lisal et al. [12]. Also for the vapour pressure good 
agreement is observed, as large uncertainties of the literature data have to be assumed if 
they have not been specified, cf. Fig. [71 bottom. 

Present 2CLJD VLE data generally show much lower statistical uncertainties than 
those of other authors, and they have the advantage that uniformly the same simulation 
method was used to produce them over a large parameter range. 

6.2 Thermodynamic consistency test 

The thermodynamic consistency of the simulation data was checked with the Clausius- 
Clapeyron equation 

dT* p%T* (1/p"* - l/p'*)' ^ ' 

The vapour pressure correlation from this work was used to evaluate the left hand side of 
Eq. (111). The right hand side of Eq. ^ was calculated from the simulation data of this 
work, the uncertainty was calculated by the error propagation law. The requirements of 
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Eq. (jlj) are fulfilled for almost all temperatures within the uncertainties of the right hand 
side of Eq. (jlj). Hence, it is concluded that the data from this work are thermodynamically 
consistent. 

6.3 Locus of the critical point 

The influence of the dipole fi*"^ and the elongation L* on the critical properties can be 
studied in Table [2] and in Figs. [1] and [21 The critical temperature T* and the critical 
density p* decrease with increasing L* at fixed dipolar momentum /i*^. They increase 
together with yU*^ at fixed elongation L*. The absolute increase AT* = T* (yU*^,L*) — 
T* (yU*^ = 0, L*) of the critical temperature with increasing yU*^ is more important for 
molecules with small L* than for those with large L*. In contrast, the relative increase 
AT* /T* [fi*"^ = 0, L*) of the critical temperature with increasing yU*^ only weakly depends 
on L*. A comparison of this finding to results from Garzon et al. who investigated 
VLE and critical data of Kihara fluids with axial pointdipole fi* and rod length L*, 
yields, that the thermodynamics of 2CLJD fluids and Kihara fluids with axial dipole are 
essentially different. Garzon et al. [36] found, that the relative increase of the critical 
temperature of Kihara fluids with axial dipole distinctly depends on the rod length L* of 
the Kihara potential. 

The critical pressure p* of the 2CLJD fluid decreases strongly with increasing elonga- 
tion L* for a fixed dipolar momentum yU*^. For a fixed elongation L*, the critical pressure 
typically shows a maximum with increasing dipolar momentum yu*^, cf. Table O This 
behaviour agrees with results from Lisal et al. E] for the 2CL JD fluid. 

As shown in Figs. [S] and HI the absolute value of the slope of the function Inp* vs. 
1/T*, i.e. the enthalpy of vaporization Ah*, increases with increasing yU*^ for a fixed 
elongation L*. It decreases with increasing elongation L* for a fixed dipolar momentum 
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/i*^. It can also be seen that the vapour pressure p* of a model fluid of given elongation 
L* decreases when the dipolar momentum /x*^ is increased. 

6.4 Deviation from principle of corresponding states 

Both the presence of a dipole and molecular anisotropy cause deviations from the simple 
principle of corresponding states [Ml [13], i.e. neither the plots ln(p*/p*) vs. T*/T* nor 
the plots T*/T* vs. p*/p* show unique curves regardless of the values of /i*^ and L*. The 
widening effect of the dipole on the density coexistence curve is shown in top of Fig. [H]for 
various 2CLJD fluids with L* = 0.6, and the displacement of the vapour pressure curves 
due to the dipole is shown in bottom of Fig. [HI These effects have already been described 
by Lupkowski and Monson [37] and Lfsal et al. [I3j for dipolar two-centre Lennard- Jones 
fluids. 

These deviations from the simple principle of corresponding states are reflected by the 
behaviour of the acentric factor uj vs. /i*^ and L*, cf. Eq. and Table [21 The more a 
fluid deviates from that principle, the higher u will be. 

A principle of corresponding states for the relative increase of the critical temperature 
has been derived by Garzon et al. [Ml [SB] for dipolar and quadrupolar Kihara fluids. 
In analogy to results in a comprehensive investigation of the 2CLJQ fluid |9|, where that 
principle of corresponding states could not be conflrmed for the 2CLJQ fluid, in this work, 
that principle could not be conflrmed for the 2CLJD fluid, either. This is a further clear 
hint on considerable differences between the thermodynamics of two-centre Lennard- Jones 
fluids and Kihara fluids. 
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7 Conclusion 



The present paper aims at the quahtative and quantitative improvement of available 
VLE data of the 2CLJD model fluid. In a systematic investigation the two parameters 
/i*^ and L* of that model fluid were varied in the ranges < /x*^ < 20 and < L* < 1.0 
respectively. A total of 38 model fluids were studied in detail, including the non-polar and 
spherical cases. The NpT+ Test Particle Method was applied for the production of VLE 
data in the temperature range of about 0.55 ■ T* < T* < 0.95 ■ T* for all model fluids. At 
the lowest temperature of systems with /i*^ < 12 the chemical potentials were calculated 
by gradual insertion. The comparison of data from this work to a hybrid 2CLJD EOS 
[HI EHl [30] reveals some shortcomings of that EOS. Critical data for all 38 systems were 
obtained from individual flts of the saturated densities with Guggenheim's equations. 

In order to obtain useful tools for adjustments of model parameters of the 2CLJD fluid 
to experimental data of real fluids, global correlations of the critical data, the saturated 
density coexistence curve, and the vapour pressure curve of the 2CLJD fluid were devel- 
oped. Data from this work agree well with results of other investigators. In most cases, 
however, data from this work have lower uncertainties, furthermore the whole data set is 
self-consistent. 

Using the correlations from this work, the influence of /i*^ and L* on the locus of 
the critical point and on the shapes of the saturated density coexistence curves and of 
the vapour pressure curves were studied. Deviations from the principle of corresponding 
states are due to the presence of polarity and anisotropy. 

The correlations resulting from this work will be used to develop molecular models for 
real dipolar fluids. 
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8 List of Symbols 



a 


interaction site counting index 


7 

b 


interaction site counting index 




coefficient of correlation function 


c 


coefficient of correfation function 


c 


constant m set of efementary functions 


77" 


TT^l 1 14- ^ ^ ^ 

Hefmfioftz energy 


T7I 


fiinction to minimize 


G 


lunction to minimize 


7 


entfiafpy 




data point counting index 


i 


efementary function counting index 


% 


particle counting index 


3 


elementary function counting index 


3 


particle counting index 


k 


Boltzmann constant 


k 


efementary function counting index 


T 

Li 


molecular elongation 





simplified notation for L* 


m 


mass of particle 


m 


simplified notation for /j, 


AT 
iV 


number of particles 


n 


number of time steps 


n 


number of Monte Carlo loops 


P 


pressure 


r 


site-site distance 


rc 


centre-centre cut-off radius 


T 


temperature 


t 


time 


u 


pair potential 


u 


internal energy 


w 


virial 


y 


linear combination of elementary functions 


z 


compressibility factor 
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a coefficient of elementary function 

P coefficient of elementary function 

7 coefficient of elementary function 

jij angle between two dipole vectors 

Ahv enthalpy of vapourisation 

ATc absolute increase of the critical temperature 

At integration time step 

S statistical uncertainty 

e Lennard- Jones energy parameter 

eg relative permittivity of dielectric continuum 

Sq permittivity of the vacuum 

9 angle of nutation 

II dipolar momentum 

II chemical potential 

^ elementary function 

p density 

a Lennard- Jones size parameter 
a standard deviation 
X elementary function 
■0 elementary function 
cu acentric factor 

Vector properties 

E electric field vector 

r position vector 

/J, dipole vector 

T torque vector 

u} orientation vector 
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Subscripts 



c 


property at critical point 


corr 


from present correlation 


u 


dipole 


estd 


estimated value 


other 


from other authors 


RF 


reaction field 


sim 


from present simulation 


a 


vapour-liquid coexistence 


2CLJ 


two-centre Lennard- Jones 


2CLJD 


two-centre Lennard- Jones plus pointdipole 



Superscripts 

* reduced 

' on bubble line 

" on dew line 

id ideal gas 

res residual property 

tot for all particles 
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Table 1: Vapour-liquid equilibrium data. Extract from simulation results for 38 model 
fluids for low, mid, and high temperatures. At the low temperatures, the second virial 
coefficient was used for the vapour phase. At the temperatures marked by f , the data 
are based on chemical potentials obtained by gradual insertion. For the remaining tem- 
peratures, the data arc based on chemical potentials obtained by Widom's method. The 
numbers in parentheses indicate the uncertainties of the last decimal digits. 





Pa 








p"* 










* 


L* = 0, n*^ = 


2.92600t 


0.00883 


(10) 


0.82998 


(20) 


0.00310 


(4) 


-26.8419 


(57) 


-0.2279 


(27) 


4.25600 


0.1469 


(14) 


0.66404 


(55) 


0.04347 


(53) 


-22.457 


(15) 


-2.459 


(33) 


5.05400 


0.4120 


(34) 


0.4925 


(31) 


0.1475 


(38) 


-17.768 


(77) 


-6.98 


(17) 


L* = 0, fi*'' = 3 


3.07290t 


0.00783 


(8) 


0.84391 


(19) 


0.00262 


(3) 


-30.0653 


(79) 


-0.2347 


(26) 


4.46960 


0.1476 


(16) 


0.67441 


(56) 


0.04122 


(46) 


-24.795 


(17) 


-2.621 


(26) 


5.30770 


0.4245 


(37) 


0.5001 


(30) 


0.1472 


(38) 


-19.507 


(77) 


-7.70 


(18) 


L* = 0, = 6 


3.33250t 


0.00599 


(8) 


0.85910 


(17) 


0.00184 


(3) 


-35.1985 


(96) 


-0.2665 


(37) 


4.84720 


0.1352 


(14) 


0.68964 


(48) 


0.03466 


(45) 


-29.143 


(16) 


-2.970 


(43) 


5.75065 


0.4053 


(35) 


0.5200 


(20) 


0.1224 


(31) 


-23.346 


(59) 


-8.11 


(20) 


L* = 0, n*^ = 9 


3.63830t 


0.00513 


(8) 


0.87156 


(21) 


0.00145 


(2) 


-40.872 


(11) 


-0.3687 


(60) 


5.29200 


0.1199 


(17) 


0.69981 


(55) 


0.02821 


(48) 


-34.055 


(20) 


-3.475 


(63) 


6.28430 


0.3909 


(35) 


0.5362 


(19) 


0.1082 


(30) 


-27.867 


(60) 


-9.37 


(22) 


L* = 0, 11*^ = 12 


4.01780 


0.00435 


(91) 


0.87717 


(26) 


0.00111 


(24) 


-46.773 


(12) 


-0.483 


(104) 


5.84400 


0.1191 


(19) 


0.69862 


(53) 


0.02602 


(39) 


-38.948 


(20) 


-4.482 


(44) 


6.93975 


0.3852 


(45) 


0.5301 


(20) 


0.0889 


(62) 


-32.038 


(67) 


-10.51 


(53) 


L* = 0, = 16 


4.53370 


0.00323 


(35) 


0.88692 


(27) 


0.00073 


(8) 


-55.163 


(19) 


-0.626 


(70) 


6.77600 


0.1352 


(28) 


0.67890 


(68) 


0.02720 


(96) 


-45.066 


(27) 


-6.84 


(29) 


7.91330 


0.4209 


(54) 


0.5062 


(31) 


0.1000 


(26) 


-37.29 


(11) 


-15.05 


(35) 


L* = 0, ^1*'^ = 20 


5.11670 


0.00305 


(48) 


0.89663 


(41) 


0.00062 


(10) 


-64.216 


(33) 


-0.94 


(15) 


7.60000 


0.1335 


(33) 


0.67349 


(71) 


0.0248 


(11) 


-52.230 


(28) 


-9.14 


(44) 


8.93080 


0.4475 


(92) 


0.4862 


(89) 


0.1029 


(41) 


-42.93 


(32) 


-19.56 


(62) 


L* = 0.2, = 


2.42000t 


0.00576 


(31) 


0.73095 


(24) 


0.00244 


(13) 


-22.5710 


(79) 


-0.1737 


(96) 


3.52000 


0.1092 


(11) 


0.58419 


(52) 


0.03850 


(56) 


-18.792 


(14) 


-2.024 


(47) 


4.18000 


0.2990 


(22) 


0.4236 


(23) 


0.1279 


(23) 


-14.595 


(52) 


-5.769 


(94) 
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Table [TJ continued. 







P'* 




p"* 










* 


L* = 0.2 




















2.53720t 


0.00506 (6) 


0.74233 


(15) 


0.00204 


(2) 


-25.1959 


(76) 


-0.1784 


(21) 


3.69040 


0.1080 (10) 


0.59307 


(52) 


0.03727 


(31) 


-20.711 


(15) 


-2.292 


(19) 


4.38240 


0.3046 (26) 


0.4289 


(35) 


0.1265 


(32) 


-15.990 


(86) 


-6.34 


(15) 


L* = 0.2, /i*^ = 6 


2.74340t 


0.00397 (5) 


0.75312 


(20) 


0.00148 


(2) 


-29.2578 


(80) 


-0.2074 


(26) 


3.99040 


0.0971 (10) 


0.60338 


(46) 


0.03053 


(47) 


-24.126 


(15) 


-2.555 


(45) 


4.73860 


0.2965 (27) 


0.4490 


(24) 


0.1115 


(30) 


-19.121 


(64) 


-6.97 


(18) 


L* = 0.2 


, = 9 


















2.98490t 


0.00315 (5) 


0.76127 


(17) 


0.00108 


(2) 


-33.6769 


(83) 


-0.2596 


(40) 


4.34160 


0.0845 (11) 


0.60947 


(45) 


0.02362 


(29) 


-27.941 


(15) 


-2.726 


(38) 


5.15570 


0.2825 (29) 


0.4630 


(18) 


0.1001 


(30) 


-22.710 


(53) 


-8.23 


(20) 


L* = 0.2, /i*^ = 12 


3.28790t 


0.00319 (5) 


0.76489 


(21) 


0.00100 


(2) 


-38.165 


(11) 


-0.3912 


(59) 


4.78240 


0.0830 (17) 


0.60746 


(47) 


0.0229 


(10) 


-31.755 


(17) 


-3.84 


(20) 


5.67910 


0.2921 (31) 


0.4581 


(16) 


0.0966 


(44) 


-25.958 


(48) 


-9.82 


(36) 


L* = 0.2 


, = 16 


















3.71250 


0.00191 (57) 


0.76778 


(17) 


0.00053 


(16) 


-44.4425 


(75) 


-0.37 


(11) 


5.44000 


0.0913 (20) 


0.59958 


(58) 


0.02294 


(83) 


-36.866 


(22) 


-5.47 


(24) 


6.41250 


0.3067 (36) 


0.4424 


(22) 


0.0966 


(26) 


-30.213 


(73) 


-12.92 


(29) 


L* = 0.4 


, = 


















1.78750t 


0.00326 (3) 


0.59571 


(11) 


0.00187 


(2) 


-17.2882 


(60) 


-0.1314 


(14) 


2.60000 


0.06524 (64) 


0.47244 


(58) 


0.03277 


(80) 


-14.213 


(14) 


-1.708 


(73) 


3.08750 


0.1892 (18) 


0.3417 


(26) 


0.1231 


(59) 


-10.955 


(54) 


-5.08 


(21) 


L* = 0.4, = 3 


1.86835t 


0.00271 (3) 


0.60327 


(11) 


0.00149 


(2) 


-19.2066 


(54) 


-0.1316 


(16) 


2.71760 


0.06292 (78) 


0.47941 


(49) 


0.0317 


(13) 


-15.655 


(13) 


-1.97 


(13) 


3.22715 


0.1830 (14) 


0.3362 


(25) 


0.1048 


(26) 


-11.808 


(58) 


-4.92 


(12) 


L* = 0.4, = 6 


2.01025t 


0.00205 (3) 


0.61042 


(12) 


0.00104 


(1) 


-22.1685 


(55) 


-0.1547 


(20) 


2.92400 


0.05441 (67) 


0.48638 


(43) 


0.02355 


(39) 


-18.139 


(13) 


-1.993 


(34) 


3.43968 


0.1632 (15) 


0.3710 


(16) 


0.0839 


(17) 


-14.631 


(42) 


-5.123 


(89) 


L* = 0.4, fi*'^ = 9 


2.17415t 


0.00156 (2) 


0.61655 


(14) 


0.00073 


(1) 


-25.4563 


(61) 


-0.1937 


(29) 


3.16240 


0.04737 (81) 


0.49242 


(42) 


0.01895 


(42) 


-20.979 


(14) 


-2.271 


(50) 


3.75535 


0.1605 (19) 


0.3697 


(16) 


0.0770 


(21) 


-16.889 


(46) 


-6.02 


(13) 
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Table [TJ continued. 









o'* 




//* 

y 






* 


j^ir res* 




L* = 0.4 


, fi*^ = 12 


















2.38315t 


0.00137 


(2) 


u.uiouo 




U.UUUOc' 


(^ \ 


98 7AnA 




n 9Rnn 

-u.zuuu 


oo j 


3.46640 


0.0455 


15) 


n 48990 


(44) 


ni 79 


('39'! 


-93 731 




-3 06 


f69l 


4.09920 


0.1590 


'28) 


n 3699 




0.0732 


(33) 


-19.401 


('471 


-7.23 




L* = 0.505, /i*^ = 


= 


















1.55650t 


0.00237 


(3) 


U.04:±U0 




U.UUIOU 


l'9'l 






n 1 1 1T? 


[10) 


2.26400 


0.05086 


'49) 




K } 




v / 


1 9 ^S"? 


(^ 1 1 




n 71 


2.68850 


0.1487 


'15) 


n 9997 




1093 


('38'l 


-9 31 


noi 


-4.42 




L* = 0.505, = 


= 3 


















1.58070t 


0.00147 


(2) 




yd) 


onoQ^ 

yJ .\j\j\JijO 




1 7 9fi9/l 


C^l 1 


u. uoy o 


n 91 


2.36500 


0.04830 


58) 






09^*^0 

yJ .\j ^OOiJ 


v / 


1 S"^^ 

±0 . OOtJ 


CI 91 


± .4:U4: 


oo J 


2.80820 


0.1479 


'18) 


0.3121 


(36) 


0.1016 


(33) 


-10.593 


C821 


-4.65 


f15l 


L* = 0.505, n*'^ = 


= 6 


















1.74520t 


0.00148 


(2) 






n nnns7 
u.uuuo ( 




1 Q RSI 9 




n 1 '?7i 

-U. iO / i 




2.53800 


0.04245 


;64) 


U.4:4:UO -L 




091 1 7 


(AOi\ 




(^ 1 1 


1 7^S 




3.01435 


0.1365 


;i2) 


n 3203 


(32) 


0.0849 




-12.472 


C751 


-4.99 


f17l 


L* = 0.505, = 


= 9 


















1.88320t 


0.00108 


(2) 


u.ooyoo 




u.uuuoy 




99 1^770 


(fu 1 


n 1 7'^i 




2.73900 


0.03763 


'74) 






n 091 1 




18^11 


(^ "^1 


9 77 


(AA 1 


3.25280 


0.1252 


'14) 


0.3321 


(16) 


0.0715 


(16) 


-14.842 


(41) 


-5.59 


(11) 


L* = 0.505, = 


= 12 


















2.06030t 


0.00096 


(1) 


0.56200 


(13) 


0.00048 


(1) 


-25.4946 


(73) 


-0.2481 


(35) 


2.99000 


0.0352 


'10) 


0.44477 


(41) 


0.0159 


(30) 


-20.985 


(13) 


-2.73 


(66) 


3.56000 


0.1256 


'19) 


0.3223 


(19) 


0.0702 


(30) 


-16.723 


(56) 


-6.79 


(20) 


L* = 0.6 




















1.40250t 


0.00197 


(2) 


0.50039 


(9) 


0.00144 


(2) 


-14.0889 


(46) 


-0.1027 


(11) 


2.04000 


0.04288 


'48) 


0.39339 


(41) 


0.02746 


(70) 


-11.408 


(10) 


-1.387 


(55) 


2.42250 


0.1286 


;i3) 


0.2643 


(34) 


0.1075 


(33) 


-8.299 


(70) 


-4.26 


(13) 


L* = 0.6, = 3 


1.40800t 


0.00106 


(2) 


0.51372 


(8) 


0.00076 


(1) 


-15.8529 


(52) 


-0.0751 


(11) 


2.12500 


0.04056 


'48) 


0.39966 


(40) 


0.02425 


(41) 


-12.580 


(11) 


-1.407 


(25) 


2.45952 


0.1047 




0.3067 


(13) 


0.0706 


(38) 


-10.124 


(29) 


-3.35 


(17) 


L* = 0.6, fi*'^ = 6 


1.56640t 


0.00117 


(2) 


0.51251 


(11) 


0.00077 


(1) 


-18.0029 


(61) 


-0.1286 


(18) 


2.27840 


0.03499 


'59) 


0.40589 


(38) 


0.01938 


(59) 


-14.577 


(11) 


-1.615 


(54) 


2.70560 


0.1181 


;i4) 


0.2987 


(32) 


0.0863 


(27) 


-11.399 


(77) 


-4.91 


(14) 
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Table [TJ continued. 







P'* 




p"* 






* 


j^ir res* 




L* = 0.6 


, = 9 


















1.68800t 


0.00082 (1) 


0.51756 


(11) 


0.00050 


(1) 


-20.6624 


(61) 


-0.1637 


(23) 


2.45500 


0.0276 (12) 


0.40916 


(39) 


0.01464 


(57) 


-16.807 


(12) 


-1.888 


(56) 


2.88515 


0.0979 (13) 


0.3160 


(11) 


0.0629 


(33) 


-13.769 


(27) 


-4.94 


(21) 


L* = 0.6, /i*^ = 12 


1.84250t 


0.00071 (1) 


0.52036 


(13) 


0.00040 


(1) 


-23.3629 


(78) 


-0.2435 


(31) 


2.68000 


0.0265 (15) 


0.40876 


(41) 


0.0123 


(13) 


-19.059 


(13) 


-2.06 


(26) 


3.15552 


0.0989 (18) 


0.3071 


(21) 


0.0632 


(37) 


-15.509 


(56) 


-6.18 


(26) 


7" * MO 

L = 0.8 


, /i ^ = 


















1.17700t 


0.00135 (2) 


0.43871 


(9) 


0.00118 


(2) 


-12.2841 


(39) 


-0.0855 


(11) 


1.71200 


0.03152 (36) 


0.34003 


(50) 


0.02282 


(22) 


-9.768 


(12) 


-1.090 


(12) 


2.03300 


0.0972 (11) 


0.2103 


(44) 


0.0969 


(32) 


-6.664 


(91) 


-3.71 


(11) 


L = 0.8 


, At = 3 


















1.17150t 


0.00065 (1) 


0.45198 


(9) 


0.00057 


(1) 


-13.8571 


(45) 


-0.05905 


(90) 


1.77200 


0.02847 (34) 


0.34700 


(35) 


0.02044 


(37) 


-10.8032 


(93) 


-1.193 


(26) 


2.10430 


0.0925 (10) 


0.2421 


(96) 


0.0952 


(36) 


-7.98 


(20) 


-4.09 


(13) 


L* = 0.8 




















1.30130t 


0.00070 (1) 


0.45104 


(10) 


0.00055 


(1) 


-15.7904 


(55) 


-0.1031 


(15) 


1.89280 


0.02375 (53) 


0.35393 


(41) 


0.01654 


(77) 


-12.575 


(11) 


-1.482 


(86) 


2.24770 


0.0846 (10) 


0.2567 


(15) 


0.0687 


(26) 


-9.692 


(38) 


-4.12 


(19) 


L* = 0.8 


, = 9 


















1.38160t 


0.00043 (1) 


0.45755 


(10) 


0.00032 


(1) 


-18.3431 


(57) 


-0.1412 


(25) 


2.02720 


0.02041 (58) 


0.36185 


(34) 


0.01305 


(52) 


-14.691 


(10) 


-1.707 


(65) 


2.40730 


0.0725 (11) 


0.2632 


(21) 


0.0517 


(19) 


-11.524 


(51) 


-4.07 


(14) 


L* = 0.8 


, 12*^ = 12 


















1.52096t 


0.00037 (1) 


0.46033 


(13) 


0.00025 


(1) 


-20.8029 


(69) 


-0.2326 


(40) 


2.19840 


0.01671 (78) 


0.36363 


(33) 


0.0084 


(11) 


-16.796 


(11) 


-1.62 


(23) 


2.61060 


0.0658 (12) 


0.2648 


(12) 


0.0442 


(44) 


-13.318 


(33) 


-4.62 


(36) 


L* = 1.0, = 


0.97900t 


0.00050 (1) 


0.41024 


(7) 


0.00052 


(1) 


-11.5972 


(43) 


-0.04220 


(67) 


1.50750 


0.02504 (30) 


0.30729 


(43) 


0.0184 


(14) 


-8.871 


(11) 


-0.931 


(89) 


1.69100 


0.05413 (62) 


0.2461 


(13) 


0.0436 


(44) 


-7.358 


(29) 


-1.93 


(18) 


L* = 1.0, fi*'^ = 3 


1.02850t 


0.00044 (1) 


0.41346 


(9) 


0.00043 


(1) 


-12.8117 


(48) 


-0.04662 


(66) 


1.50750 


0.01711 (27) 


0.32581 


(31) 


0.01378 


(28) 


-10.1499 


(90) 


-0.848 


(17) 


1.77650 


0.05688 (68) 


0.25218 


(96) 


0.0495 


(52) 


-8.122 


(23) 


-2.39 


(24) 



29 



Table [TJ continued. 





Pi 




p"* 




* 




+ 


L* = 1.0, fi*'^ = 6 


1.08900t 


0.00026 (1) 


0.42225 (8) 


0.00024 (1) 


-15.0377 


(51) 


-0.0568 


(11) 


1.62000 


0.01403 (48) 


0.32977 (34) 


0.01014 (65) 


-11.761 


(10) 


-0.921 


(74) 


1.90000 


0.05034 (73) 


0.2575 (10) 


0.0464 (16) 


-9.533 


(27) 


-2.98 


(10) 


L* = 1.0, = 9 


1.17700t 


0.00018 (1) 


0.42930 (9) 


0.00015 (1) 


-17.5873 


(90) 


-0.0976 


(45) 


1.71200 


0.01020 (51) 


0.34335 (28) 


0.0103 (18) 


-14.056 


(10) 


-2.24 


(56) 


2.03300 


0.04132 (86) 


0.26785 (87) 


0.0294 (28) 


-11.437 


(24) 


-2.77 


(22) 


L* = 1.0 


, fi*^ = 12 














1.30740t 


0.00014 (1) 


0.43197 (11) 


0.00011 (1) 


-20.166 


(12) 


-0.1711 


(40) 


1.90160 


0.01033 (65) 


0.34090 (33) 


0.00665 (27) 


-15.996 


(12) 


-1.484 


(15) 


2.28190 


0.0501 (10) 


0.2482 (23) 


0.0471 (31) 


-12.506 


(62) 


-5.09 


(32) 
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Table 2: Critical data, the critical compressibility factor, and the acentric factor of 38 
2CLJD model fluids. 





L* 









0.2 


0.4 


0.505 


0.6 


0.8 


1.0 






n 


5.236 


4.313 


3.163 


2.735 


2.454 


2.049 


1.762 


rp* 

c 






0.3113 


0.2710 


0.2251 


0.2032 


0.1850 


0.1577 


0.1153 


0* 






0.4736 


0.3670 


0.2116 


0.1616 


0.1353 


0.0995 


0.0794 








0.2878 


0.3106 


0.2972 


0.2908 


0.2980 


0.3079 


0.3101 


Zc 

c 






-0.0542 


-0.0245 


0.0308 


0.0514 


0.0589 


0.0755 


0.1382 


UJ 






5.475 


4.521 


3.311 


2.871 


2.558 


2.127 


1.876 


rp* 
C 






0.3197 


0.2808 


0.2237 


0.2028 


0.1868 


0.1611 


0.1449 


P* 






0.5157 


0.3619 


0.2133 


0.1711 


0.1378 


0.0932 


0.0856 








0.2946 


0.2851 


0.2880 


0.2939 


0.2884 


0.2720 


0.3149 


Zc 






-0.0191 


0055 


0528 


0.0720 


0958 

\J * KJ \J \J KJ 


0.1475 


0.1797 


CO 






5.990 


4.917 


3.591 


3.100 


2.774 


2.316 


2.010 


c 






0.3169 


0.2781 


0.2202 


0.2017 


0.1839 


0.1576 


0.1428 


pI 

rc 






0.5179 


0.3624 


0.2105 


0.1677 


0.1342 


0.0894 


0.0824 


v* 

rc 






0.2728 


0.2650 


0.2662 


0.2682 


0.2631 


0.2449 


0.2871 


Zc 






0.0214 


0.0471 


0.0976 


0.1179 


0.1429 


0.1987 


0.2398 


UJ 




9 


6.585 


5.382 


3.909 


3.377 


3.040 


2.520 


2.178 


C 






0.3124 


0.2700 


0.2164 


0.1959 


0.1790 


0.1510 


0.1382 


Pi 






0.5098 


0.3594 


0.2078 


0.1646 


0.1304 


0.0839 


0.0765 


P*c 






0.2478 


0.2473 


0.2457 


0.2488 


0.2396 


0.2205 


0.2542 


Zc 






0.0687 


0.0938 


0.1457 


0.1669 


0.1936 


0.2565 


0.3114 


UJ 




12 


7.289 


5.945 


4.282 


3.676 


3.273 


2.713 


2.370 


rp^ 

C 






0.3033 


0.2647 


0.2122 


0.1921 


0.1765 


0.1519 


0.1376 


Pi 






0.5067 


0.3619 


0.2095 


0.1651 


0.1291 


0.0797 


0.0716 


Pi 






0.2292 


0.2300 


0.2306 


0.2338 


0.2235 


0.1934 


0.2196 


z. 






0.1159 


0.1396 


0.1921 


0.2145 


0.2435 


0.3163 


0.3881 


UJ 
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Table [21 continued. 





L* 









0.2 


0.4 


0.505 


0.6 


0.8 


1.0 




r' 


16 


8.249 


6.654 












T* 

C 






0.2961 


0.2565 












Pr 






0.5202 


0.3795 












Pr 

Jr c 






0.2130 


0.2224 












Zc 






0.1732 


0.1945 












UJ 




20 


9.164 














C 






0.2884 














Pi 






0.5567 














Pi 






0.2106 














Zc 






0.2220 














UJ 
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and L* = 0.3292 ( ); [I3] 2CLJD with p*^ = 12 and L* = 0.22 ( ); 
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Appendix A 

Technical details 

The two-centre Lennard- Jones plus axial pointdipole (2CLJD) pair potential writes as 

wherein M2clj is the Lennard- Jones part 

2 2 

M2CLj(ni, uJi,u}j,L) = ^^Ae 

a=l 6=1 

and md is the dipolar part, as given by Gray and Gubbins 

1 /x^ 

uuirij, UJi, ujj, i/) = ^ (cos7ij - 3cos6'iCos6'j) . (5) 

4vreo \rij\ 

Herein is the centre-centre distance vector of two molecules i and j, Tab is one of the 
four Lennard- Jones site-site distances; a counts the two sites of molecule i, b counts those 
of molecule j. The vectors LOi and ujj represent the orientations of the two molecules i and 
j. 6i and 6j are the angles between the dipole vectors /Xj and of the molecules i and j 
respectively and their centre-centre distance vector Vij. The cosine of the angle between 
the dipole vectors /Xj and /x^- is calculated as cos7jj = (/Xj ■ / (l^tj ■ \ The number 
of parameters related to the dipole is reduced to one, namely the dipolar momentum /i, 
as its position and orientation within the molecule are fixed and as it is reduced by the 
large distance approximation to a pointdipole whose interaction is described by Eq. ([5]). 
The Lennard- Jones parameters a and e represent size and energy respectively. 

Beyond a certain elongation L, for small intermolecular distances \rij\ the M2Cljd pair 
potential might diverge to infinity, as the positive Lennard- Jones term M2clj cannot out- 
weigh the divergence to infinity of the dipolar term md, that occurs for some relative 
orientations of the molecules i and j. This divergence leads to infinite Boltzmann factors, 
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i.e. non-existence of the configurational integral. During phase space samphng by molec- 
ular simulation within the pressure range in question, this artefact of the 2CL JD potential 
causes no problem as intermolecular centre-centre distances are very improbable to fall 
below critical values. However, the calculation of the chemical potential by test particle 
insertion often runs into critical intermolecular centre-centre distances. To avoid compu- 
tational problems in such cases, the dipole site was shielded by a hard sphere of diameter 
0.4 ■ a. This shielding of the dipole site is in analogy to the suggestion of MoUer and Fis- 
cher [23] to shield the quadrupole site in two-centre Lennard- Jones plus pointquadrupole 
fluids, where similar problems occur. This hard sphere was not active during configura- 
tion generation. For reasons of consistency, this shielding by a hard sphere with diameter 
0.4 ■ a was applied to all fiuids studied here. 

In order to achieve a monotonous transition from L > to L = the spherical 
fluids ICLJ, where /i^ = 0, and ICLJD (Stockmayer fluid), where L > were treated 
as two-centre LJ fluids with L = 0. This leads to site superposition that is not present 
when ICLJ and ICLJD fluids are represented by a simple LJ site. Therefore the reduced 
temperatures, reduced pressures, reduced enthalpies and reduced dipolar momenta here 
are fourfold of the corresponding values in the one LJ site case. Densities, of course, are 
not concerned. 

For all simulations, the centre-centre cut-off radius was set to 5.0 ■ cr. Outside 
the cut-off sphere the fluid was assumed to have no preferential relative orientations of 
the molecules, i.e., in the calculation of the LJ long range corrections for the potential 
internal energy, the virial, and the chemical potential, orientational averaging was done 
with equally weighted relative orientations as proposed by Lustig [10]. Long distance 
corrections for the dipolar part of the potential model were calculated with the reaction 
fleld method [^ W2\ . The reaction fleld is derived from the polarization of the dielectric 
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continuum supposed to surround the cut-off sphere of molecule i and is calculated as 

For sufficiently large systems, i.e. N > 500, the sensitivity of simulation results to the 
value of the relative permittivity eg is negligible, cf. Saager et al. [12]. Therefore, the 
relative permittivity eg is usually set to infinity in simulations of dipolar fluids, cf. [I3l [HI 
HH]. For reasons of consistency, in this work, the relative permittivity eg was uniformly 
set to infinity. 

The interaction of the reaction field -Erf,* with the dipole /x, of molecule i contributes 
to the potential energy of this molecule 

The dipolar potential internal energy per particle of the system of 2CLJD particles is 
then 

N-l N 



4 = 1 J=i + 1 

rij<rc 



, 1 2 (eg - 1) 1 



1 es - 1 ] 



Aireo 2eg + 1 r 

Knowing m^*, the dipolar contribution to the virial is easily calculated by 

w'^' = -3 ■ <\ 

The homogeneous reaction field -Erf.i in each cut-off sphere exerts no force on the dipole 
in the centre of the cut-off sphere, however, its contribution to the torque on molecule 

i is 



TRF,i — Hi X -ErF.i, 

which has to be added when molecular dynamics simulations are considered. 
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Appendix B 

Correlation of critical properties 

The development of correlations in the present work was done in analogy to a previous 
work on the 2CLJQ fluid |9]. The functions T*{fi*^,L*) and pl{fi*'^,L*) were assumed to 
be linear combinations of elementary functions one of which is a constant c, the others 
depend either on /i*^, i.e. ipilfi*"^), or on L*, i.e. C,i{L*), or on both, i.e. Xt if^*"^ ^ L*) . 
The number of elementary functions was restricted to up to two for both the /i*^- and 
the L*-dependence and to up to three mixed terms. With y representing any of the 
aforementioned functions, the linear combination writes as 

<2 <2 <3 

y L*) = c + |] a. ■ + J2f^r ^ALl + E 7. ■ Xk{^i*\ L*). (6) 

i=l j=l k=l 

Usual non-weighted least squares fits of these linear combinations to the critical data 
given in Table [2] yielded the coefficients in Table El which also contains the elementary 
functions. It should be mentioned, that the elementary function C,i (L*) from [9j could be 
reused here, whereas new elementary functions ipi and Xt were selected due to different 
macroscopic thermodynamics of dipolar and quadrupolar fluids. 

Correlation of saturated densities and vapour pressures 

The saturated density - temperature correlations are based on Eqs. ([H) and ([2]). The 
dependence on /i*^ and L* was ascribed to T*, p*, and to the coefficients Ci to C3. The 
coefficient functions to C3 (/i*^, L*) were linear combinations of elementary 

functions in the sense of Eq. 1^. The correlations T* (/i*^, L*) and p* (/i*^, L*) were used 
in the fit of the functions p'*{p*^, L*, T*) and p"*{p*^, L*, T*) to data from simulation. The 
set of elementary functions used in the coefficient functions L*) to C3 (/i*^, L*) were 

different from those used for the 2CL JQ fluid [9] . 
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For the vapour pressure - temperature correlation the polynomial ansatz 

lnp^(/i ,L,T) = ci{fi + — + — , (7) 

was applied in analogy to [9J. The coefficients Ci(/i*^,L*) to c^{jj*'^,L*) were linear com- 
binations of the elementary functions in the sense of Eq. These elementary functions 
were different from those used for the 2CLJQ fluid [9]. 

The fit of the functions p'*{fj.*'^, L* ,T*), p"*{fx*^, L*,T*), and Inp^p*^, L*,T*) to data 
from simulation was performed by uncertainty- weighted least squares minimization of the 
functions F and G defined as 



E 



f*2 



I 

= min, 



2 ! 

= mm. 



The resulting coefficients and the elementary functions are given in Table El 
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Table 3: Elementary functions and their coefficients for the correlations T* {jj,*"^, L*), pi (/x*^, L*),Ci (/x*^, L*) to C3 (/x*^, L*), 
and ci (/x*^, L*) to C3 (/x*^, L*). The notation is simplified: i is L*, m is /x*^. 



r*(TO,£) 




c 


1 


0.1454013-101 


c 


1 


0.3157828 




TO/(88 + m)2 


0.1363894-103 




m/(88 + 777)2 


0.9871123- lOi 


toV(88 + m)3 


0.2020243 • 10^ 


7n2/(88 + m)3 


-0.1461751 ■ 103 




1/(0.1 + ^2) 


0.3269772 




£2/(0.11+^2) 


-0.1475616 


1/(0.1 + ^5) 


0.4910240-10-1 




-0.4152214- 10-1 




m/((88 + m)2. (0.1+^2)) 


0.4239005 - 102 


Xt 


777/(88 + m)2 .£2/(0.11 +£2) 


-0.1010584- 102 


mV((88 + m)3 • (0.1 + £2^) 


0.6724083 - 10^ 


77l2/(88 + 77l)3 -£2/(0.11 +£2) 


0.4105884- 102 


mV((88 + m)3.(0.1+^5)) 


0.7913876 • 102 


m2/(88 + m)3. £5/(0.11 +£5) 


0.5299302 ■ I02 






c 


1 


0.2951644 


c 


1 


0.6484789 - lO-i 


ipi 


m2/(70 + m)2 


-0.6339151 




m'^/{70 + 7n)2 


0.7301440 


771^/(70 + 771)3 


0.3182745 - 10^ 


771^/(70 + 77l)3 


-0.8780100-101 






-0.2359527 


^^ 


£2.e^ 


-0.6551324 


£3 


0.5466755 


£3 


0.1810641 • lOi 




m2/(70 + m)2 .f7(£ + 0.4) 


0.1449170- 10^ 




7772/(70 + 777)2 .£2 .e*^ 


0.4808117- lOi 


m^/{70 + mf -i^/ii + OAf 


-0.1955388 


777-7(70 + 777)3 .£2 .e^ 


0.1937455-101 


toV(70 + to)3 ■£8/(£ + 0.4) 


-0.5849357- 10^ 


7^2/(70 + 777)2 - £3 


-0.1320822- 102 


C3 




c'i 


(m,£) 


c 


1 


-0.7258204-10-2 


c 


1 


-0.5486341 • 10-2 




m2/(70 + m)2 


-0.6215183 




7772/(70 + 777)2 


0.1223952- lOi 


771^/(70 + 771)3 


0.4708560 - 10^ 


7773/(70 + 777 )3 


0.1350701 - IQi 






0.4316296 


^^ 


£2/(0.4 + £)2 


0.2479957 


£3 


-0.1166922-101 


£7(0.4 + £) 


-0.1684560 




7n/(88 + m)2.fV(0.4 + ^)2 


-0.2215021 ■ 102 


Xt 


m/(88 + m)2. £2/(0.4 + £)2 


-0.1956742 • 102 


mV(88 + m)3 .^2/(0.4 + ^)2 


-0.7803181-102 


7772/(88 + 777)3 .£2/(0.4 + £)2 


-0.2289032- 103 


rn/(88 + ?B)2 .£8/(0.4 + £) 


-0.5735507 


m2/(88 + m)3 -£7(0.4 + £) 


0.7221121 - 103 



Table [3l continued. 





ci (m,£) 


c 


1 


0.2574709-10-1 


c 


1 


0.4411718- 10* 




m/(88 + m)2 


-0.2940407-101 




m/(88 + TO)2 


0.4575129- 103 


mV(88 + m)3 


-0.1008706-103 


m2/(88 + TO)3 


0.2469929 ■ 10'' 






-0.9426323-10-1 


^^ 


£2/(0.4 + £)2 


-0.2016356-101 




0.1108324 


£V(0.4 + £) 


0.4346103 




m/(88 + TO)2 -^7(0.4 + ^)2 


0.2543158-102 


x% 


m2/(88 + m)3 -£2/(0.4 + £)2 


0.9787962 - 103 


mV(88 + m)3 .^2/(0.4 + 


0.5987224-102 


m2/(88 + m)3 -£8/(0.4 + £) 


0.2467171 • lO'' 


m7(88 + m)3 .^7(0.4 + 


-0.6511462-103 








C3 (m, £) 


c 


1 


-0.2686327-102 


c 








m/(88 + TO)2 


-0.3428826-10-* 




m/(88 + TO)2 


-0.5264689 • 103 


mV(88 + m)3 


-0.8720808- 10^ 


m2/(88 + TO)3 


0.6782756 • W 




£2/(0.75 + 


0.1275315-103 




t 


0.1812550 


£3/(0.75 + ^)3 


-0.1393077-103 








m/(88 + m)2 .£2/(0.75 + ^)2 


0.7248855 - lO"* 








m2/(88 + m)3 . £2/(0.75 + ^)2 


0.5715498- 10^ 






m2/(88 + mf ■ £3/(0-75 + £)3 


-0.7433962 - 10^ 







